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Abstract. We present new computations of approximately length-minimizing polygons with fixed 
thickness. These curves model the centerlines of "tight" knotted tubes with minimal length and fixed 
circular cross-section. Our curves approximately minimize the ropelength (or quotient of length and 
thickness) for polygons in their knot types. While previous authors have minimized ropelength for 
polygons using simulated annealing, the new idea in our code is to minimize length over the set of 
polygons of thickness at least one using a version of constrained gradient descent. 

We rewrite the problem in terms of minimizing the length of the polygon subject to an infinite 
family of differentiable constraint functions. We prove that the polyhedral cone of variations of a 
polygon of thickness one which do not decrease thickness to first order is finitely generated, and 
give an explicit set of generators. Using this cone we give a first-order minimization procedure and 
a Karush-Kuhn-Tucker criterion for polygonal ropelength criticality. 

Our main numerical contribution is a set of 379 almost-critical knots and links, including all 
prime knots with ten and fewer crossings and all prime links with nine and fewer crossings. For 
links, these are the first pubhshed ropelength figures, and for knots they improve on existing figures. 
We give new maps of the self-contacts of these knots and links, and discover some highly symmetric 
tight knots with particularly simple looking self-contact maps. 



1. Introduction 



Overview. Knots tied in rope are flexible machines which organize tensions and contact forces 
to bind tightly and resist unravelling. As a technology, knots have proved remarkably effective. 
For this reason there is a vast body of knowledge about their practical uses. Yet in many ways, 
the design of these machines remains mysterious. As early as 1987 Maddocks and Keller were 
able to study different types of hitches and predict their holding power by an analysis of their 
equilibrium shapes [34|. But these shapes were rather simple, and there was no way to infer the 
structures of more complicated knots from these examples. It was obvious that what was needed 
was data, and by the end of the century a series of numerical experiments in knot-tightening were 
underway [[44ll4n[32l[52ll . This paper describes a new computational approach to knot-tightening 
which yields improved numerical results (a preliminary report on some of our findings appeared in 
the conference proceedings [13J). To build our method, we derive some new results in the theory 
of ropelength for polygonal knots. 
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Defining the problem. Given any space curve 7, we can define the thickness Thi(7) of 7 to 
be the supremal e for which any point in an e-neighborhood of 7 has a unique nearest neighbor 
on the curveQ Any curve with nonzero thickness is C^'^ (that is, is with a Lipschitz first 
derivative) [l22l[TTI . Given this, it has been shown that 

Proposition 1.1 dlSBl ). If •y is a curve, then the thickness Thi(7) is given by the supremal 
radius of all embedded tubes formed by taking the union of disks of uniform radius centered on 
7(5) in the planes normal to 7'(s). 



This idea of thickness was first proposed by Krotenheerdt and Veit in 1976 [|30l [3T1l and was 
rediscovered in the 1990's by Nabutovsky [37J and Buck and Orloff [|6J. The thickness can be used 
to define a scale-invariant quantity called ropelength: 

Definition 1.2. The ropelength of a curve 7 is defined by 

where Len(7) is the length of 7. The minimal ropelength of a knot or link type L, Rop(L), is the 
minimal ropelength of all curves in that knot or link type. 



The knot tightening problem is to find and describe the minimal ropelength curves in a given 
knot type. It is known that such curves exist, but their exact shapes are currently the subject 
of active mathematical research (c.f. [|25l l26l [TTII '). Once found (or computed to sufficient accu- 
racy), these configurations have been used to predict the relative speed of DNA knots under gel 
electrophoresis [|28l . the pitch of double helical DNA [|36l . the average values of different spatial 
measurements of random knots Il20ll . and the breaking points of knots [|42l . They also provide a 
model for the structure of a class of subatomic particles known as glueballs flU. 



Another form of the problem. Let 7 : 5^ — t- now be a parametrized curve, and define the 
self-distance function d : 5^ x 5"'^ — t- M of 7 by d{s, t) := \\'j{s) — 7(t) ||. As usual, let k(s) denote 
the curvature of 7. We then define the set dcsd(7) of doubly-critical self-distances to be the set of 
critical points of d with s ^ t. Taking the partial derivatives of d, we see that (s, t) E dcsd(7) if 
and only if 

(7(s) -7(t),7'(5)) = Oand (7(5) - 7(t), 7 '(t)) = 0. 

A key idea in [|33ll is that for any r < Thi(7), the surface of the tube of radius r around 7 has no 
self-intersections and is smooth. But when r = Thi(7), the tube is pinched or has a tangential 
self-intersection. This leads to an alternate characterization of thickness: 

Theorem 1.3 ( [[33l ). The thickness of'y is the minimum of 

.1 . d{s,t) 
mm and mm . 

s k{s) (s,t)edcsd{7) 2 

Figure [T] shows curves where the first and second of these terms control the thickness, 
^ederer referred to this number as the reach of 7 ll22l . 
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Figure 1 . The thickness of a smooth curve 7 is controlled by curvature (as in the 
left picture), and the length of chords in dcsd(7) (as in the right picture). 



Since length and thickness scale together, minimizing ropelength is the same as minimizing 
length over the set of curves with thickness at least one. Since thickness is a min-function, the 
condition Thi(7) > 1 can be viewed as an infinite family of inequality constraints on 7. These 
constraints are active at places where the tube around 7 forms kinks (where 1/k is in control of the 
minimum in Theorem 1.3) or has self-contacts (where the self-distance '^{«'*)/2 is in control of the 
minimum). 



Numerical approaches to the knot-tightening problem. Previous authors have defined dis- 
cretized versions of thickness for polygons or spline curves and viewed the problem as one of 
minimizing the nonsmooth quotient of length and thickness. The advantage of this approach is 
that it is a very simple and robust way to obtain approximately ropelength-minimizing curves. The 
disadvantage is that it is very difficult to take advantage of the fact that thickness (as given in 
Thm|1.3|) is a min-function. 



Our approach is to define a discrete version of thickness as a min-function and think of the 
problem as one of minimizing a differentiable function Len(V) subject to a family of differentiable 
constraints Thip(V) > 1. While our approach will not quite fit into the standard framework of 
constrained optimization (our family of constraints is infinite), we will be able to define a version 
of constrained gradient descent which minimizes polygonal ropelength effectively. 



Theoretical framework. For an equilateral space polygon V we first prove that our function 
Thip(V) can be written as a min over a fixed compact family of differential functions. From here 
we use Clark's theorem to show that Thip has a one-sided derivative in the direction of any variation 
W of V. For a polygon with Thip(V) = 1 we use these derivatives to define a cone of infinitesimal 
variations I(V) which do not decrease Thip to first order and the dual cone of "resolvable" vari- 
ations R(V). Our next main theorem is that R(V) is a finitely generated polyhedral cone whose 
generators are the gradients of the lengths of certain chords of the polygon (called struts) and of a 
function of certain turning angles of the polygon (called kinks). We give explicit formulae for these 
gradients in terms of the vertex positions. We then compute the gradient of Len(V) and define the 
constrained gradient of length to be the projection of Len(V) onto the polyhedral cone I(V). At 
this point we give the expected result that a polygon is critical for polygonal ropelength if and 
only if the constrained gradient of length is zero. Equivalently, a polygon is critical for polygonal 
ropelength if there is a set of positive Lagrange multipliers on the struts and kinks which combine 
to equal the negative of the length gradient. The theory section ends with a discussion of how to 
compute the constrained gradient numerically. 
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Numerical methods. Sections [3] and |4] describe the design of our polygonal ropelength minimiz- 
ing software. Our algorithm essentially consists of computing the constrained gradient of length 
and taking small steps in this direction until the constrained gradient is sufficiently small. However, 
the details of the process are not quite so simple. Since the constraint functions are nonlinear, even 
steps that are in the direction of the constrained gradient violate some constraints to second order. 
Further, newly active constraints are discovered throughout the run as previously distant sections 
of tube come into contact with one another. As a result, we must choose stepsizes carefully and 
correct errors periodically. It is also important to run efficiently, as the size of our problem (about 
one thousand variables and a similar number of active constraints) is fairly large. We have solved 
these technical and engineering problems and used our software to minimize all prime knots with 
ten or fewer crossings and all prime links with nine or fewer crossings, for a total of 379 different 
knot and link types. We intend to address the ropelength of composite knots and links in a future 
publication. 

New ropelength bounds. We check our figures against previous computations of the minimum 
ropelength of knots and links and against some of the few known theoretical results for the lengths 
of tight links. Our results improve on all previously published computational results except for 
the trefoil knot. For example, we improve the best known upper bound for the ropelength of 
the well-studied figure-eight knot 4i by 0.06 to 42.0887 (as compared to the bound of [|T4ll ) and 
improve the best known upper bound for the ropelength of the 92o knot by 8.12% to 80.2219 
(compared to the bound of [4Tj). To get a sense of the difference between the configurations 
produced by our method and the configurations produced by the simulated annealer of [47] we 
show both configurations in Figure |2] For links, our figures are the first computational results 
to appear in print, but compare well to known theoretical results. For example, the upper bound 
provided by our computation of the Borromean rings link 62 is 58.0070 — within 0.0017% of the 
exact value around 58.0060 suggested by [10], while our computation of the tight shape of the 
"simple chain" link is 41.7086588 — within 0.02% of the correct value of of 67r + 2 [11,1. 

We also compared our results to those of Gilbert [|24|. which are unpublished but available on 
Bar-Natan's Knot Atlas wiki. Gilbert provides Fourier cofficients and instructions for reconstruct- 
ing the vertices of his configurations from this data. We followed his instructions, but our software 
did not verify his claimed ropelength number^ According to our measurement of the ropelength 
of Gilbert's configurations, our knots are tighter in all cases but 2f by an average of 3.619%, with 
some outliers, such as our 9^7 link, which is 71% shorter. If we compare our results to Gilbert's 
claimed ropelengths, our knots and links are tighter in 316 cases and less tight in 36. Overall, our 
knots and links are (on average) 1.075% tighter than the bounds claimed by Gilbert with our 928 
Unk about 4% shorter than Gilbert's claim. 

Self-contact maps. Two sets of authors (von der Mosel et al. [|50l and Cantarella et al. ifTOl ) have 
given versions of a ropelength criticality criterion for knots without kinks which state roughly that 
a knot 7 is ropelength-critical when the elastic force given by the gradient of the length of the curve 
is balanced by a system of Lagrange multipliers on the self-contacts of the tube around 7. The latter 

^Our measurement of curvature by MinRad is sensitive to edgelength and seems to come out much larger than his 
ropelengths would indicate. This is probably a discretization effect, and it is certainly very possible that the Fourier 
knots defined by Gilbert's data have ropelengths corresponding to Gilbert's claimed numbers. 
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Figure 2. These two images of the 92o knot show the tightest configurations ob- 
tained by our algorithm (left) and by the TOROS algorithm described in [|47l (right). 
It is clear that our algorithm performs better once there are many self-contacts in the 
knot. In fact, the ropelength of the left-hand configuration is bounded by 80.2219, 
while the configuration on the right has ropelength bounded by 87.31. 

authors used their condition to derive a ropelength critical configuration of the Borromean rings 
and a surprising ropelength critical configuration of a clasp formed by two tubes stretched across 
each other. 

In each of these examples, the most difficult part of the result was the deduction of the structure 
of the set of self-contacts for the tight configuration. Since these contact maps are very sensitive to 
small perturbations of the centerline, it has been difficult to resolve them using previous numerical 
method^ These contacts and the system of Lagrange multipliers on them are explicitly computed 
by our algorithm, allowing us to give medium-quality contact maps for a large number of knots 
and links. The contact maps offer some support for the hypothesis that a relatively small number 
of structures may reappear often in tight knots and links. 

Previous work. This is not the first time gradient-like methods have been attempted for the knot- 
tightening problem. Our work has been inspired by Piotr Pieranski's SONO algorithm ll4T]| . which 
follows a version of the length gradient, but does not include an explicit resolution of this vector 
against the active constraints. Our thinking is also informed by John Sullivan's "energy-ropelength 
method" [52J, which optimizes thickness instead of length, estimating the maximum diameter 
of a uniform embedded tube around the core curve by an average of the radii of embedded 
cross-sectional disks and minimizing the resulting smooth functional using the conjugate-gradient 
implementation in Brakke's evolver Q. 



The notable exception to this rule has been the "biarc" spline-annealing method of lfT4l . which has produced 
well-resolved contact maps for the 3i and 4i knots. 



2. A DISCRETIZATION FOR THE ROPELENGTH PROBLEM 



Polygonal thickness. Consider a closed space polygon V with vertices vi, . . . ,vy and edges 
ei, . . . , ey. We will think of V as the vector (f i, . . . , vy) in {E?Y = , and assume that all 
subscripts on vertices and edges are taken mod V . The unit tangent vector Tj to each edge of a 
polygon is well-defined on the interior of the edge. At the vertex f , joining edges ej_i and e,, 
there are two tangent vectors Tj_i and Tj. The curvature of V at Vi is usually thought of as a delta 
function whose mass is given by the turning angle 6*, from Tj_i to Tj. We will use a somewhat 
different definition of curvature for polygons: 

Definition 2.1. The minimum radius of curvature (or MinRad) of V at Vi is given by the radius of 
the unique circle that is tangent to the two edges meeting at Vi and that touches the midpoint of the 
shorter one. 

Rawdon has shown [|44| that if 6i is the turning angle of V at f j, then we can give MinRad (fj) 
(and define MinRad^ (fj)) by the expressions: 

<" '"^WW^ ^ [WFY ^ m,„{M,„Rad-(..,, M,„Rad-(„.,}^ 

It is clear that while MinRad t>j is not neccesarily a differentiable function, the two functions 
MinRad^ vi are differentiable when they are defined. The motivation for this definition is that we 
can round off all the comers of V by splicing in these circle arcs, generating a C^'^ curve with radii 
of curvature equal to the MinRad(f j). We could have defined Thip(V) to be the thickness of this 
curve. It turns out, however, that there is no closed form computation for that number (though it 
can be computed approximately, as we will see in Section [5]). 

We now define a set corresponding to dcsd for polygons: 

Definition 2.2. Let dcsd(V) be the set of (p, q) on V with p q which are local minima of the 
self-distance function on V. 

There are several possible cases for (p, q) in dcsd(V), since the polygon might have a vertex at 
one or both of the endpoints of the chord. These are shown in Figure [3j 




Figure 3. We see three types of local minima of the self-distance function on a 
space polygon V in the three-dimensional drawings above. From left to right, these 
are an edge-edge pair, a vertex-edge pair, and a vertex-vertex pair. 
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We can then define Rawdon's polygonal thickness: 

Definition 2.3. The polygonal thickness Thip(V) of a space polygon V without self- intersections 
is given by the minimum of 

d{p, q) 



Thip(V) := min < minMinRad(fi), min 

i (p,(?)edcsd(V) 

We have carefully constructed this definition so that when polygons Vn with increasing numbers 
of edges are inscribed in a space curve 7 under some mild geometric hypotheses, Thip(V„) — )■ 
Thi(7) [gailllllTl. 



The problem with Thip. Definition 2.3 allows us to define the set of polygons with Thip(V) > 1 
as the polygons obeying a family of constraints in the form MinRad(f j) > 1 and d{p, q) > 2 for 
{p, q) G dcsd(V). This is almost the standard form for constrained optimization problems: 

(2) ^miri^ /(V) subject to giiy) > 0, 

where / and the gi are differentiable. The problem is that the set of constraint functions d{p, q) for 
(p, q) G dcsd(V) depends on the polygon. We will need a common set of constraint functions for 
all polygons in a neighborhood of a solution. 

Constraint thickness. To solve this problem, we will define a new thickness measure for poly- 
gons called constraint thickness which is given in the form above. We will then prove that for 
equilateral polygons, the new constraint thickness defines the same set of polygons as the old 
polygonal thickness. 

We first define a subset of the pairs of points on a polygon 

Definition 2.4. For a given positive r and £, let 6{t, i) be the turning angle of a pair of edges of 
length i with MinRad = r. We set 

VB(r,£) = |(p,g) G V X V : vb(p,g) > ^ 

where vb(p, q) is the smaller number of vertices between points p and q (counting p and/or q if 
they are vertices) if they are on the same connected component of V and 00 otherwise. 

We note that an easy computation shows that 9{t, £) = 2 arctan(y2T). We can now define our 
new thickness measure 

Definition 2.5. The (r,^)-constraint thickness CThi(r, i, V) of a polygon V is given by 

/ON n-ru-f e • / • MinRad(t;i) . d{p,q) 

(3) L im(r, Vj = mm < mm , mm 

I T (p,g)eVB{r,Q 2 



We note that V need not be equilateral or have edgelength i to define the constraint thickness 
to defined the constraint thickness of V. We can view r as the "stiffness" of the rope (c.f. the 
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definition of A-thickness in fH and 0), as it provides a lower bound on the radius of curvature 
of a tube of unit radius. Though our theory (and our code) should work for any r > 1, we have 
not experimented with values for r other than 1 and so will write the (1, £) -constraint thickness 
CThi(l,£,V) as CThi(£,V). 

We can now prove that CThi(£, V) is an equivalent thickness to Thip for equilateral polygons of 
edgelength i. 

Theorem 2.6. If V is an equilateral polygon of edgelength i, Thip(V) > 1 <^=^ CThi(£, V) > 1. 

To prove the theorem we will need a lemma (c.f. Lemma 13 of B6l ): 

Lemma 2.7. IfV is an equilateral polygon of edgelength i and MinRad > r, then dcsd(V) C 
VB(r,£). 

Proof. The proof has two parts — in the first, we show that the shorter arc between any (p, q) ^ 
VB(r, tj has total curvature t less than vr, while in the second we will show that any pair joined by 
such an arc cannot be in dcsd(V). So suppose that t > tc. We will prove that (p, q) G VB(r, i). 

Since MinRad(V) > r, we know that each turning angle of V is less than 9{t,£). If the total 
curvature of the arc joining p and q is at least vr, then vh{p, q) ■ 6{t, i) > vr, so 

(4) "''^"•"I'^WJ, 
and (p, q) G VB(r, i), proving the claim. 

Now suppose that (p, q) G dcsd( V) . We claim that the total curvature t of each arc joining p and 
q is at least vr, and hence that (p, q) G VB(r, t). Suppose not. The arc of V joining p and q together 
with the chord from pio q form a closed space polygon V'. The total curvature of this polygon is 
equal to t plus the turning angles at p and q. By Fenchel's Theorem [fT9l . that total curvature is 
at least 2tx. So the angle at p and the angle at q must sum to more than tt. Thus either the angle 
at p or the angle at q must exceed But in that case, we could reduce d{p, q) to first order by 
moving p ox q along an edge from the arc which connects p and g, contradicting our assumption 
that G dcsd(V). □ 



We are now ready to prove Theorem 2.6 



Proof. Suppose that CThi(^, V) > 1. This implies that min^ MinRad(f j) > 1 by the definition of 



CThi. Lemma 2.7 tells us that dcsd(V) C VB(1, tj, so we know that 



(5) min d(p,q)> min d(p,q). 

(p,q)edcsd{V) (p,g)GVB(l,^) 

Together, these facts imply that Thip(V) > 1, proving one direction of the theorem. 



Suppose that Thip(V) > 1. As above, this means that min, MinRad(fi) > 1, so Lemma 2.7 



applies and ([5]) holds. If the minimum in the right-hand side of (|5]) is achieved on the interior of 




Figure 4. The key step in the proof of Theorem 2.6 is the proof that points and 
g' on an arc V are at least distance 2 apart. This arc has equal edgelengths £, each 
turning angle equal to 6{l,t) := 2arctan(y2) and n := [7r/6'(l,£)] edges. We 
see above that these conditions imply that V has an inscribed circle of unit radius. 
Further, the marked point q' must have a larger y-coordinate than the top of the 
circle, providing the required lower bound on the distance from p' to q' . 



VB(1, £), then it is a local minimum of d{p, q) where "p ^ q and so is in dcsd(V). In this case, ([5]) 
is an equality and CThi(£, V) > 1, completing the proof. 

We are left with the case where the minimum of d{p,q) over VB(1,£) is realized by some 
(p, q) on the boundary of VB(1, £). We claim that a!{p.«')/2 > 1. This will complete the proof that 
CThi(£,V) > 1. 

By definition, (p, g) is on the boundary of VB(1,£) only if vb(p, g) = \'^/e{i/)\. And since 
vb(p, q) is constant on the interiors of edges, one of p and q (without loss of generality, q) must be 
a vertex. Since each turning angle of the arc of V between p and q is bounded by ^^(1, tj, Schur's 
theorem [15J implies that d{p, q) is bounded below by the distance between the endpoints of p' , q' 
of a planar polygonal arc V with the same edgelengths and each turning angle equal to 6'(1, £). We 
depict the situation in Figure |4} 

We know that V has n = vh{p, q) edges and total curvature [n — l)6{l,tj. Since n = vh(p, q) = 
, we have 

(6) n-l< <n so {n - l)e{l, i) < n < ne{l, i). 

0[1, q 

Thus if we add an edge to V at q' with turning angle 9{l,i)to form an arc V^, the total curvature of 
V is less than vr while the total curvature of is at least tt. These facts imply that if the first edge 
of V lies along the x-axis, the point q' has the largest y coordinate on V^. But our turning angle 
and edgelength conditions imply that P+ has an inscribed circle of unit radius, so the y-coordinate 
of q' is at least two. This implies that d{p', q') > 2, completing the proof. □ 



These proofs imply an obvious corollary which will be useful in practice: 
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Corollary 2.8. 7/"dcsd(V) C VB(r, t) and the distance between any two vertices on the boundary 
o/ VB is strictly greater than Thip(V), then CThi = Thip /or polygons in a neighborhood ofV 
(regardless of whether or not V is equilateral with edgelength i). 



Proof. The argument is the same as that of Theorem|2.6[ using the hypotheses instead of Lemma|2.7 



and the argument about turning angles. □ 



Struts and Kinks. In our definition of Thip, we saw that pairs of points in dcsd and vertices 
with minimum MinRad were in control of thickness. We now want to develop similar sets of 
"controlling" pairs of points and vertices for CThi. This will require a bit of care. 

Given any two line segments ei and 62 in space, a calculation reveals that the minimum distance 
between them is attained at a single point unless ei and 62 are parallel. In that case, the minimum is 
attained at an interval of corresponding pairs (as in Figure |5]). The endpoints of these intervals are 
self-distances measured from an endpoint of one segment to a point on the other. Following this 
line of argument we see that for any space polygon the local minima of the self-distance function 
d{p, q) are isolated unless there are pairs of parallel edges, in which case there may be families of 
local minima as above. Using these observations we define 

Definition 2.9. The strut set Strut (V) is the set of pairs (p, q) in VB(1, i) with '^(p.9)/2 = 1 and 
either 

• (p, q) is an isolated local minimum of d{p, q), or 

• (p, q) is an endpoint of a family of local minima of d{p, q). 

In the second case, (p, q) must be a vertex-edge pair joining two parallel edges of V. 

We note that Strut (V) is a finite subset of dcsd(V) (which may be infinite if two edges are 
parallel). It is much easier to define 

Definition 2.10. The kink set Kink(V) is the set of vertices Vi and signs ± with MinRad^ Vi = 1. 

The strut and kink sets are both empty if CThi(£, V) > 1. 

Polygon space and variations of CThi. We now want to describe the space of variations of a 
polygon which preserve or increase CThi to first order. Given a polygon V G M^^ we can define a 
variation of V by any W = {wi, . . . , wy) E M^^. This variation generates a family of polygons 

(7) Vt = V + tW = {vi+ twi, ...,vv + twv). 

We now want to prove that CThi(£, V) has a one-sided derivative as we vary V according to any 
variation W and to give a finite procedure for computing that variation. This will require some 
setup. 
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S Vj 



Vi-i p Vi = q 



Figure 5 . When the edges Cj and ej are parallel, many chords realize the minimum 
distance between the segments. In this case, we show that the minimum derivative 
of distance between any of these pairs occurs at one end or the other. We name the 
endpoints of this family of chords p and q on and r and s on ej. One of each of 
these pairs must be an endpoint — in this case it is g = Vi and r = vj-i that are 
endpoints. 



Proposition 2.11. Suppose that CThi(£, V) = 1. Then viewing every pair of points {p,q) on V 
and every MinRad^ Vi as functions oft, the forward time derivative below exists and satisfies 



(8) Dw CThi{i,V) 



dt- 



CThi(Vi) 



mm 



t=o 
d 



mm (MinRad=^ (t) 



mm 



d d{p{t),q{t)) 



t=o 



Strut(V) dt+ 



t=0. 



Proof We begin by ignoring any MinRad Vi functions which are not defined (which happens when 
Vi^i, Vi and f j+i are colinear). Since CThi(£, V) is equal to 1, the MinRad of these vertices will not 
affect CThi(V + tW) for small enough t. The function CThi is then the minimum of a set of dif- 
ferentiable functions MinRad^ Vi and '^iP'i)/2 indexed by the (compact) disjoint union of compact 
sets {vi, ±} U ■ • ■ U {vv, ±} U VB(1, 1) (where we assume that any Vi with MinRad undefined 
are missing). Clark's theorem for min-functions [.16.1 tells us immediately that the derivative in ([8]) 
exists. 



However, Clark's theorem tells us that 



Dw CThi(£, V) = min 



mm 

{V^,±) 



d 
dt + 



(MinRad^ fi)(t), min 



. MinRad =t Vi=l 

MinRad=^ 



t=o 



{p,g)GVB{l,^) dt+ 
d(v,q)/2=l 



d{p{t),q{t)) 



t=0 



The first set {(z, ±) | MinRad Vi = 1} is the kink set, which matches ([8]). But if a pair of edges 
in V are parallel and at distance 2 from one another, then Strut (V) is only a subset of {(p, q) E 
VB(1, £) I dip,qy2 = 1}. We must prove that 



(9) 



mm 



d dipit),qit)) 



(p,g)eVB(l,^) dt+ 

d{p,ij)/2 = l 



mm 



d dipit),qit)) 



t=o 



(p,g)eStrut(V) d)f:+ 



t=0 



For any pair of parallel edges with distance 2, we may assume that the situation is as in Figure |5| 
We label points p, q, r and s as in the Figure, and parametrize the line segments between p and q 
and between r and shy r] G [0,1]. The pairs with i] = and = 1 are in the strut set of V, but the 
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pairs given by all other values of i] are not. To prove (|9]) we must find 

d 

min — — \\t]p + (1 — 'qjq — r]r — [1 — 'i])s\\ , 
»7e[o,i] dt^ 

and show that it is attained at r/ = or r/ = 1. If we view p, q, r, and s as functions of time, then 
for any given r], the time derivative of the corresponding length is given by 

-{r]P + (1 — v)q — w ~ ~ v)^^ vp' + (1 ~ v)q' ~ w' ~ (1 ~ 'n)'^')^ 

where we have used the fact that <^'y^i^^i)j2 = 1. Regrouping, we can rewrite this as 

^{r]{p - r) + (1 - r]){q - s),r]{p - r)' + (1 - r]){q - s)'), 

and using the fact that p — r = q — s at time 0, we can again rewrite this as 

ri{p — r,p — r') + — f]){q — s, q — s'). 
Now as 7] varies between and 1, we note that the r] derivative of the above quantity is 

{p — r,p' — r') — {q — s,q' — s') . 

In particular, this derivative is nonzero for all r] G [0, 1] unless {p — r,p' — r') = {q — s,q' — s'), in 
which case it vanishes identically. This means that the minimum value of this expression is always 
realized when r] = or rj = 1. This completes the proof. □ 



We can use Proposition 2.1 1 to define two sets of variations that will be of particular interest to 
us. The first set consists of variations that are tangent to the boundary or pointing into the interior 
of the set of polygons CThi(£, V) > 1. We will allow our polygons to move in these directions. 

Definition 2.12. Suppose we have a polygon V and a variation W of V. If CThi(£, V) = 1, we 
say W is an infinitesimal motion of V if the forward directional derivative 

(10) DH/CThi(^,V) > 0. 

If CThi(£, V) > 1, we call every variation W an infinitesimal motion. The set of all infinitesimal 
motions of V is denoted /(V). 



The following Corollary follows directly from Proposition 2.1 1 



Corollary 2.13. The set /(V) is the dual cone of the set — V'^(P''?)/2 /or (p, g) G Strut(V) and 
-V MinRad^ Vi for {v^ ±) G Kink(V). 



Proof. We need only recall that the dual cone to a set of vectors A is the set of vectors X for 
which (X, W) < for all W e A. Since the directional derivatives of d{p,q)/2 and MinRad"^ Vi in 
the direction X are the dot products of X with —^d.ip^i)/2 and —V MinRad^ Vi, X is in the dual 
cone if and only if all these directional derivatives are nonnegative. But by the Proposition, this 
implies that Dx CThi(^, V) is nonnegative as well. □ 



The second set of variations of interest will be the normal cone of the boundary of the set of 
polygons with CThi(£, V) > 1. We will forbid our polygons from moving in these directions. 
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/(V) 



Figure 6. The infinitesimal motions /(V) and the resolvable motions RiV) of 
V form dual convex cones. Hence, although these are not orthogonal subspaces 
of M^^, a similar decomposition property holds true: any vector W may be written 
uniquely as a sum of a vector W[ G /(V) and a vector Wr E -R(V). 

Definition 2.14. The convex cone of resolvable motions -R(V) of V is the cone generated by the 

set -V^(P'i)/2 for {p, q) G Strut(V) and -V MinRad^ Vi for [vi, ±) G Kink(V). R{V) is the set 
of vectors R G M?^ which can be expressed in the form 

(11) R= Yl -A-V^^+ Yl -A-VMinRadt;,. 

(p,<j)eStrut(V) VjeKink{V) 

Here the indices i and j just number the elements of the strut and kink sets. The constants and 
are nonnegative numbers, as suggested by the notation. 

It is a standard fact from optimization theory that -R(V) = I{V)^, since for any set of vectors 
{v} the double dual is the cone generated by {v}. 

Theory of constrained optimization. Given a function /(V) on the space of polygons M'^^, we 
can compute the negative gradient —V/, which is a variation vector in M^^. We are now interested 
in understanding how this gradient is modified by the constraint CThi(£, V) > 1. This thickness 
constraint models the effect of an embedded tube around the polygon: it allows some motions of 
V and blocks others. 

Definition 2.15. The constrained gradient (—V/)/ of — / is the closest vector in /(V) to — V/(V). 

We now recall that any convex cone and its dual cone provide a kind of orthogonal decomposi- 
tion of their ambient vector space, as shown in Figure [6} 

Proposition 2.16 dlSTI . Thm. 2.8.7). Any vector W G R^^ may be uniquely written 

(12) W = Wr + Wi, where {Wr, IV/ ) = 0, 
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Wr G -R(V) is the closest resolvable motion to W, and Wj G /(V) is the closest infinitesimal 
motion to W. 



We note that this Proposition shows that the constrained gradient of — / is well-defined. Further, 
it is easy to show that the constrained gradient is the direction of steepest descent for / within /(V). 
This makes us guess that the constrained gradient should vanish at a critical point for minimizing /. 
To prove it, we define critical points more carefully 

Definition 2.17. We say that V is thickness-critical fijr minimizing f if either: 



• Dwf = 0, or 

• CThi(£, V) = 1 and for any W with DwfiV) < 0, we have Dw CThi(£, V) < 0. 



In the first case, we are at an unconstrained critical point of the objective function /. In the 
second, we are at a constrained critical point where motion in the direction of the negative gradient 
of / is blocked by active constraints. We then have a version of the Kuhn-Tucker theorem (restated 
in our language from the original form in filOil ). which gives a verifiable condition for thickness- 
criticality. 

Theorem 2.18. The polygon V is thickness-critical for minimizing f <^==^ —V/ is in -R(V) <^==^ 
the constrained gradient (—V/)/ vanishes. 



Proof. It suffices to show that the first two statements are equivalent, since the second and third are 



clearly equivalent by Proposition 2.16 



If —V/ is not in -R(V), then Farkas' theorem implies that there exists some W with {W, V/) = 
Dwf < and {W,R) < for all R G R{V) (L40J, p. 118). Using the definition of RiV) and 
Proposition 2.1 1 , this implies Dw CThi(£, V) > 0. Thus V is not thickness-critical for minimiz- 
ing /. 



If —V/ is in -R(V) we will prove that V is thickness-critical for minimizing /. We first observe 
that the dual cone of —V/ contains the dual cone R^{V). Now suppose we have some W with 
Dwf < 0. Then {W, -V/) >0,soW ^ (-V/)+ and in particular W ^ i?+(V). But this means 
that {W, R) > for some R G RiV), so Dw CThi(£, V) < 0. Hence V is thickness-critical for 
minimizing /. □ 



We can give a natural interpretation of this Theorem in mathematical and physical terms by 
considering the condition —V/ G RiV). By definition, this means that 

(13) -V/+ Yl \'V^^+ Yl A,2VMinRadt;, = 0. 

(p,g)GStrut(V) i.jeKink(V) 

Mathematically, the Af and A^ are Lagrange multipliers. If we think of the thickness constraint as 
an embedded tube around V, we can interpret these scalars as magnitudes of compression forces 
transmitted by tube contacts (for struts) and angles where the polygon resists further bending (for 
kinks). 
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In general, we cannot expect every local minimum of a constrained function to be a constrained 



critical point in the sense of Definition 2.17 If the set of polygons defined by CThi(£, V) had an 
outward-pointing cusp we might reach a point where some W with Dwf < had D]y CThi = 0. 
For example, the constrained system 

minimize f{x, y) = — x, subject to g{x, y) = min{x^ — y,y} > 

has this property at the local minimum (0, 0) for W = (1, 0). The problem here is simply that 
Dw9 < for all W. This does not happen for thickness-constrained polygons, but we will need 
another idea to prove it: 

Definition 2.19. We say that V is constraint-qualified (in the sense of Mangasarian and Fro- 
movitz 11351) if there exists some W so that Dw CThi > 0. 

It is then standard to show 

Proposition 2.20 ( ifTOl '). Any constraint-qualified local minimum of f is a thickness-critical point 
for minimizing f. 

In our case, scaling V provides the desired motion, so we have 

Corollary 2.21. If the polygon V is a local minimum for f, then it is a thickness -critical point for 
minimizing f. 

We make a final note that in general, our criticality theory works equally well for CThi and 
Thip (even for polygons V which are not equilateral), as long as they obey the hypotheses of 



Corollary 2.8 This is true in practice in all of our numerically computed configurations. 

3. Bridging Theory and Computation 



Overview of the algorithm. We have now derived enough theory to describe our algorithm in 
general terms. We wish to minimize the function Len(V) subject to the constraint CThi(^, V) > 1. 
We will do so by computing the constrained gradient (— Len V)/ and stepping in this direction. 
These steps will reduce Len(V) while keeping V close to the set CThi(£, V) > 1 (since the con- 
straints are nonconvex, we cannot stay entirely inside this set). When (— Len V)/ vanishes, the 
algorithm will terminate. By Theorem [2.18 if the constrained gradient was exactly zero, the re- 



sulting configuration would be a thickness-critical point for minimizing length. We note that our 
algorithm will attempt to maintain an approximately equilateral polygon V but it is not required 



to: constant edgelength i is not a hypothesis of Theorem 2.18 Our only caveat is that we must 



remember that CThi(V) may not be equal to Thip(V) if the final configuration fails to obey the 



hypotheses of Corollary 2.8 We also note that there is nothing special about choosing Len(V) as 
the function to minimize — both our theory and our code would work just as well for any other 
function. 



Computing the constrained gradient. To implement this algorithm, we must be able to compute 
the constrained gradient (—V/)/. This is a standard problem in linear algebra. By definition, if 
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—V/ is written as (— V/)/j + (—V/)/ using Proposition 2.16 the constrained gradient is equal 
to (—V/)/. We can compute that by computing (— V/)ij, which is easy to do since we know the 
generators of the cone -R(V). 

Definition 3.1. If CThi(£, V) = 1, the rigidity matrix A of V is the matrix whose columns are the 
gradients -V'^(P''?)/2 for {p, q) G Strut(V) and -V MinRad"^ Vi for {vi, ±) G Kink(V). 

We can construct the rigidity matrix by finding the members of Strut (V) and Kink(V). It 
follows from the definition that -R(V) is the image of the positive orthant under the matrix A. By 



Proposition 2.16 (— V/)/j is the closest vector in that image to —V/. So if we solve the non- 



negative least- squares (NNLS) problem 

(14) mm||/lA- (-V/)||, 

then {—Vf)R = AA and (—V/)/ = —V/ — AA. This least-squares problem is a special kind of 
quadratic programming problem which has been well-studied in numerical linear algebra (see [4J). 
In our case, the problem is much easier because A is extremely sparse — the gradients of the <^(.P'<i)/2 
involve no more than 4 vertices (and so 12 variables), while the gradients of the MinRad^ involve 
only 3 vertices (and 9 variables). So each column of A, which is typically 1000 or more entries 
long, contains at most 12 nonzero entries. 

The gradient of Length. We can now compute (—V Len)/ if we can compute —V Len, build the 



rigidity matrix A from the strut and kink sets, and solve the NNLS problem in ( [14] ). We will take 
these problems in order. 

Length is a differentiable function of polygons V G M^^, whose gradient is given by a straight- 
forward calculation: 

Proposition 3.2. The gradient of length of a polygon Vn is given by the collection ofn vectors 

(15) VLen(V)fc = T^ rr + j. rr. 

\\Vk-l-Vk\\ ||^^fc+l-ffc|| 



The gradient of '^(P'9)/2. Given a pair of points {p, q) on V, the gradient of the distance between 
them is a set of four vectors located at the endpoints of the edges on which p and q lie. These 
vectors are given by a calculation: 

Proposition 3.3. Suppose that {p,q) G Strut(V). Ifp = avi + (1 — and q = (3vj + (1 — 

j+i then 

= {a{p -q),{l-a){p-q), - p), {1 - P){q - p)} . 

where these three vectors are applied to vi, Vi^i, vj and Vj+i in order 

The gradient of MinRad^. As we noted above, the MinRad^ are differentiable where they 
are defined. We now compute the gradient on MinRad'*', noting that the gradient of MinRad^ is 
similar. 

16 



Proposition 3.4. Given a vertex i on Vn with finite MinRad^(f »), we let n denote the oriented 
normal vector to the plane defined by f fj, f j+i and define the scalar constant 

^ \\Vi+i-Vi\\ 



2cos0 - 2 
and the vector constants 



2 



Then if we write the gradient o/MinRad"*" as a triple of vectors located at Vi-i, Vi, and fj+i we 
have 

VMmRad+(fi) = {W, -W - X - V, X + V}. 



Proof. The proof is a lengthy calculation. We want to compute the gradient of MinRad^(t>j) = 
, where 9 is the turning angle at vertex Vi. We start with a change of variables. Let 
Vi and B = Vi^i — Vi. We can rewrite MinRad^ in terms of these variables and 



2 tan(e/2) 

A = Vi-i 
compute its gradient as follows: 



(16) 

Now 

(17) 



151 



2 tan(^'/2) 2tan(y2) 



B 
\B\ 



1 



151 



d 



tan(^/2) 



I l+cos i 



1 + COS 6* 



tan2(e/2) d9 

and tan^(^'/2) 



1 — cos 6 
1 + cos 6* ■ 



So we can rewrite (fT6l) as 
V 



151 



1 



2 tan(V2) 2tan(V2) 



0, 



B 

\B\ 



\B\ 



2 - 2 cos 



= {0,V) + KV9. 



Keeping track of the sign of the exterior angle, we see that if n is the oriented unit normal to the 
plane containing A and B, we have 



A X n n X B 



B 



\B\ 



so 



2 tan(^'/2) 



{W,X + V). 



Using the definition of A and B to change back to the original variables completes the proof. □ 



The function MinRad(f j) provides a discrete analog to the radius of curvature for the polygonal 
curve V at fj. Since this is a numerical computation of a second derivative, we expect the function 
to be quite sensitive to small changes in the positions of the vertices of V. This sensitivity will 
limit the accuracy of our computations, so we record an estimate of the norm of the gradient of 

MinRad+(t;,). 

Corollary 3.5. IfV is an equilateral polygon with edgelength i and MinRadfj = 1 then 

VMinRad^-yJ > — . 
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Proof. Consider 

_ 11^7^+1 - Vi\\ \\{vi-i - Vj) X n|| 
|2cos^ — 2| — 

Since the polygon is equilateral, and n is a unit vector normal to — Vi, this is just = 
y|2cos6i-2|. If MinRad = 1, then (squaring MinRad and using both half-angle formulae for tan- 
gent) we see that || = |2+2cose|/£2 Since W appears alone in the formula for V MinRad^, this 
is a lower bound for the norm of the entire gradient. □ 



4. Program design 



Issues of scale. The design and implementation of our algorithm ridge runner were shaped 
by the scale of the knot-minimizing problems we intended to solve and the amount of computer 
power we had on hand to solve them. To inform the discussion that follows, we will now take a 
moment to consider the dimensions of our problems. In a typical run, we started by minimizing 
the length of a low-resolution version of our knot or link with 2 vertices per unit of ropelength (80 
to 150 vertices). Once that configuration was minimized, a medium resolution run at 4 vertices 
per unit of ropelength was performed. A final run followed at 8 vertices per unit ropelength. 
Most of the runtime was spent during the final run, which took 20 — 40 CPU hours on a desktop 
computer. During the final run, the average edgelength £ for our curves was approximately 0.061, 
which meant that there were 658 edges. The average size of the strut set was 919 pairs of points, 
while the average size of the kink set was 19 vertices. The rigidity matrix was then on average 
a 938 X 1974 matrix which was 99.4% sparse (no more than 11199 of its 1851612 entries were 
nonzero). A typical run contained several hundred thousand steps. 



The algorithm. Our method is based loosely on the method of constrained gradient descent. 
The basic idea is to generate a series of polygons Vj which converge to a limit polygon which is 
thickness-critical for minimizing a function /(V) by taking a series of steps in the form 

(18) Vfc+i = Vfc -t- «(— V/)/, where a is chosen by a search algorithm. 

When CThi(^, V) > 1, this is just the method of steepest descent, since (—V/)/ = —V/. When 
CThi(£, V) = 1, these steps are tangent to the boundary of CThi(£, V) > 1 and in principle 
decrease CThi by no more than 0{a^). In some circumstances, such as when two sections of tube 
touch for the first time, we can decrease CThi by 0{a) (which is much larger, since a « 1). We 
control this error by searching for an a which keeps CThi(£, Vk + a(— V/)/) within acceptable 
bounds. When CThi(£, V^) becomes too small, we correct the accumulated error using a Newton's 
method-type solver. The code terminates when we the constrained gradient is small enough to 
convince us that we are near a point which is thickness-critical for minimizing /. This procedure 
is summarized in Algorithm [1] 

In the rest of this section, we will comment on each of these steps in turn. 
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Algorithm 1: The outline of the ridgerunner algorithm. 



input : A polygon Vo and an error bound error bound MaxErr. 
output: A sequence of positions Vk with CThi(£, Vfc) > 1 - MaxErr. 

1 repeat 



2 Compute -V/ = -V Len(Vfc) + -V Eq(Vfc); 

3 Find Strut (V) and Kink(V) and construct the rigidity matrix A; 

4 Compute constrained gradient (—V/)/.; 

5 Search for a so that Vk + a(— V/)/ minimizes ropelength and is computationally 
acceptable and set Vfc+i = + a(— V/)/; 

6 if CThi(£, Vfc+i) < 1 - MaxErr then 

7 Correct CThi(£, V^+i) by Newton's method; 

8 end 



9 until ll(^'^/)-rll/!|-v/|| is sufficiently small; 



Step 2. Equilateral polygons, CThi and Thip. We have only proved that CThi > 1 <^=^ 
Thip > 1 for equilateral polygons. It is therefore important that our Vk remain at least approxi- 
mately equilateral during a run. We enforce this constraint by defining a penalty function Eq(V) 
which is minimized when Vk is equilateral and minimizing the sum Len(V) + Eq(V). This is 
quite effective (a typical run recorded an average error in edgelength of about 0.385%) in practice. 
We note that while CThi and Thip might not be equal for nonequilateral polygons, we avoid any 
problems that might result by performing all of our final ropelength calculations with respect to 
the original Thip thickness. 



Step 3. Finding Strut (V) and Kink(V). In principle, the strut and kink sets could be found by 
direct inspection of all pairs of edges and all vertices of V. But since there are usually 10^ such 
pairs, this naive method consumes too much runtime. So to find the strut and kink sets, we used 
the clustering code oct rope of Ashton and Cantarella described in [H]. This was fast enough that 
over 30 seconds of a typicaj^run about 10% of runtime was spent finding Strut (V) and Kink(V). 
The algorithm in oct rope does not take advantage of the fact that it is called successively on data 
which vary little between calls, so a much faster customized strut-finding code could be written 
into ridgerunner. However, these figures show that this project would have little impact on 
overall performance. 



Step 4. Finding the constrained gradient. Once we have Strut (V) and Kink(V) we can use the 



gradient formulae given in Propositions 3.3 and 3.4 to construct the rigidity matrix A. We must 
then solve the sparse non-negative least squares (SNNLS) problem minA>o||^A — (— V/)||, which 



we recall as Equation 14 on page 16 



We use the freely available tsnnls library of Cantarella, Piatek, and Rawdon [[T2]| . which 
is an implementation of the block-pivoting algorithm of Portugal, Judice and Vicente ^43\ . The 
PJV algorithm solves a sequence of unconstrained least-squares problems to find a partition of 



"^a 400 edge 5.1 knot with about 600 struts 
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the variables of A into complementary sets F and G representing variables which will be nonzero 



and zero in the solution to ( [T4| ). It is very important to take advantage of the sparsity of A in 
order to solve these (rather large) problems in an acceptable amount of time, as this step makes 
the dominant contribution to our overall runtime in most cases. To this end, tsnnls solves the 
least-squares problem Ax = bhy solving the "normal equations" Ax = A^h. Since A is 
symmetric, we can solve this system using a Cholesky factorization. This is done very quickly 
using the multifrontal supernodal sparse Cholesky code TAUCS of Toledo etal. fl53l . 

We have sacrificed some accuracy in favor of speed, since the condition number of A is the 
square of the condition number of A. A standard "rule of thumb" in such situations is that the 
error in the solution is on the order of machine epsilon (10~^^) multiplied by condition number. 
To verify that this was small in practice, we used the rcond function in LAPACK to estimate the 
condition number of the rigidity matrices of all of our final configurations. The average condition 
number was on the order of 10"^ with none being worse than 8 x 10^. Thus we expect to have an 
average error on the order of 10^^ and a worst-case error of 10^® in our final computations of the 
constrained gradient. 



It is also worth noting that the TAUCS code will fail if the rigidity matrix is singular, which 
will occur when there is more than one way to balance gradient force. This is expected for very 
complicated knots, but seems to be rare among knots in our dataset. A more advanced version of 
tsnnls would calculate a minimum-norm solution to the least-squares problem in this case. 



Step 5. Choosing a stepsize. When CThi(V) > 1 our code sets a small maximum stepsize 
of 10^^ and proceeds by Euler integratiorj^ Once CThi(V) = 1, thickness typically decreases 
by a small amount on each step. We choose a by a line search algorithm, finding the minimum 
ropelength of configurations in the given direction using Brent's method with a relatively low 
precision. 

However, we do not always accept the ropelength-minimizing a. Instead, we apply a collection 
of ad hoc conditions which we describe as a being "computationally acceptable". These include an 
upper bound on stepsize of 10^^, a lower bound of 10^^, and the requirement that the linear algebra 
solver of Step 4 can compute a new direction —V// at the new location. These are motivated by 
several practical considerations. If the stepsize is permitted to be too large, loose configurations 
will often form large kinked regions before the tube contacts itself. Kinks reduce stepsizes by 
orders of magnitude- in practice, this means that such a run takes an unacceptably long time to 
converge. If the stepsize is permitted to be too small, the solver can stall just before discovering a 
new self-contact. In these cases it has proved better to take the risk of a slight increase in ropelength 
in order to improve the strut set. Finally, even when the stepsize is less than 10~^, if an arc of the 
knot suddenly contacts another arc, introducing too many new struts into the rigidity matrix, the 
matrix can become numerically singular, defeating the tsnnls solver of Step 4. Thus, we must 
look ahead and make sure the next position will be acceptable to tsnnls before locking in a 
stepsize. 

^We could improve the accuracy and speed of this portion of the computation by using a smarter ODE solving 
method. But these steps have no linear algebra involved, so they are already orders of magnitude faster than the ones 
to come. In practice, this portion of the run consumes < 1% of the total runtime. 
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Step 7. Error correction. When the error bound MaxErr = 10~^ is reached, we use Newton's 
method to return Vk to a configuration with larger thickness. For any given variation of V we can 
estimate the change in the <^(p.'?)/2 for {p, q) G Strut(V) and in MinRad^ Vi for (f j, ±) e Kink(V) 
by A^W , where A is the rigidity matrix we have already computed. 

We use this observation in a straightforward way. We construct a vector C of desired corrections 
which is equal to (1 - MaxErr/2) - a!{p,g)/2 for (j9, q) G Strut (V) and (1 - MaxErr/2) - MinRad=^ Vi for 
(fj, ±) G Kink(V). Having done so, we find a minimum-norm solution to A^W = C. We then 
step according to W, using a search algorithm to decide the stepsize, rebuild the rigidity matrix in 
case we have changed the strut or kink set in the correction step, and iterate. 

We note that we do not attempt to correct all of the error in CThi(V) during this procedure. If 
we did so, we would risk losing struts and kinks when we rebuild the rigidity matrix. In that case, 
the next Newton step, ignoring those pairs or vertices, might rediscover them as struts and kinks. 
In principle, this cycling behavior could delay or prevent convergence of the Newton procedure, 
as noted by Fletcher [23J. Our method does not eliminate this possibility entirely (in the current 
version of the code, we have observed occasional failures of the Newton solver) but in practice the 
Newton solver almost always converges in only a few iterations. 

The main problem with the Newton solver is that it is slow for large problems. The matrix 
is mapping from a high-dimensional space of variations to a relatively low-dimensional space 
of struts and kinks, so it has a large kernel. Hence the matrix AA^ is not positive definite, and 
so we cannot solve A^W = C using the method of normal equations and the fast Cholesky 
decomposition of TAUCS. Instead, we must use the older Isqr code of Paige and Saunders (39] 
to find a minimum-norm solution to the problem. This can be very slow. For instance, in a 640 
edge trefoil with 975 struts and 10 kinks, correction steps consumed anywhere between 3 and 25 
seconds of runtime. Normal steps completed in less than a second. We always have the option of 
sidestepping Newton correction by simply scaling the knot (as in Pieranski's SONO algorithm). 
This preserves ropelength but destroys the strut set completely, requiring us to rebuild the strut set 
during subsequent steps. Our experience has been that this can improve performance during the 
middle stages of a run, when a fairly large number of struts and kinks have formed but the knot is 
still far from tight, but it is better to use Newton correction in the final stages of a run when one is 
trying to adjust a converged strut set to improve the final results. 

At the moment, the speed of Isqr controls the overall performance of our code. We hope to 
find an improved error-correction procedure in future versions of the software. 



Modified versions of the algorithm. We have also modified our algorithm to handle some 
special cases, such as open curves with fixed endpoints or endpoints constrained to lie in planes. 
In these cases, the gradients of the endpoint constraints are added to the rigidity matrix and the 
gradient of length is resolved against them in Step 4. In addition, a specialized error-correction 
algorithm enforces the constraints after each step to prevent numerical error from causing the 
endpoints to drift away from their positions over time. The general Newton's method algorithm 
for error-correction is also modified in these cases to take endpoint constraints into account. 
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In addition, we have found that curves whose final tight positions have long segments with 
no struts or kinks as well as tightly curved regions with many struts and kinks often take a very 
large number of steps to tighten completely. Sections of the curve with no struts or kinks simply 
minimize length with no constraints and must therefore end up as straight lines. But as they 
approach this position, the gradient of length approaches zero, while regions where the gradient 
of length is balanced by struts and kinks have comparatively large length gradients. Since the 
step size is controlled by the tightly curved regions, it may take a very long time for the strut and 
kink-free regions to finish straightening. We have had some success in these cases with a modified 
version of our algorithm which detects sections of curve with no struts or kinks and scales up the 
length gradient on those portions of the curve alone. 



5. Results of Computations 

We now present the main results of our computations. To summarize, we have significantly 
extended the range and quality of existing computations of tight knots and links. The new data 
support some interesting conjectures about the geometric structure of these configurations. 

Validation of ridgerunner computations. To verify that the system works, we checked the 
results of ridgerunner against some theoretical results. The results of the comparison appear 
in Table [T| As we can see from the Table, the relative error in these ropelength computations is as 
small as 0.0017%. 

The paper (W\ also gives an explicit strut set for the Borromean rings. To compare the numeri- 
cally computed strut set to the theoretical one, we plot them together in Figure|7J The Figure shows 
that the numerically computed strut set is quite close to the actual one. Figure [8] shows a similar 
comparison between theoretical results and a ridgerunner computation for the strut set of the 
"simple clasp" formed by two strands looped over one another. The theoretical results in [fTOl for 
this clasp assume that the curvature of the clasp is not bounded, so we compare with the results of 
a run of our software which did not enforce curvature constraints. 

Computing polygonal ropelength minimizers for many knots and links. We minimized polyg- 
onal ropelength for all prime knots of 10 and fewer crossing and all prime links of 9 and fewer 
crossings (a total of 379 knot and link types) at resolutions of at least 8 vertices per unit of rope- 
length (several hundred vertices in total). For a few knots and links of special interest, we computed 
high resolution runs with 16, 32, or 74 vertices per unit ropelength. The largest runs in our dataset 
contain about 2400 vertices. 

The computations were performed on clusters at the University of St. Thomas, the University of 
Georgia, and the ACCRE cluster at Vanderbilt University. We began our computations with an ini- 
tial low-resolution (200 vertices or fewer) polygon, which we ran until the residual was sufficiently 
low. We then increased resolution by a minrad-preserving version of spline interpolation and mini- 
mized again from the resulting new starting configurations. Our inital goal was a residual less than 
0.01, which we achieved for 375 of the 379 knots and links in our data set. We were able to reach 
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29.01 ^ 



26.10. 



23.20. 



19.34. 



^ 20.30. 



17.40. 



14.50. 





Figure 7. The diagonal above is labeled with arclength values along the three 
components of the Borromean rings link, which is numbered 62 in Rolfsen's table. 
Every pair G Strut(V) is represented by a dark green square centered on 

{s,t). As we see from the top plot, no tube around a component of the link is in 
contact with itself (so the three triangles near the diagonal are empty). But each of 
the components makes contact with the other two, as shown by the boxes plotted in 
the rectangles forming the remainder of the plot. We can see that the contacts break 
up naturally into "lantern- shaped" structures. In the bottom plot, we compare one 
"lantern" to the self-contact set predicted by [10], which is represented by a black 
line. 
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Link name 


Clasp 


Hopf link (21) 


2?#2? 


Borromean rings (63) 


Vertices 


332 


216 


384 


930 


Ropp bound 


4.2841 


25.1406 


41.7131 


58.0192 


Rop bound 


4.2837 


25.1334 


41.7086588 


58.0070 


Smooth length 


4.2629[fT0l 


Stt^ 


127r + 4|fTll 


58.0060[|TOl 


Relative error 


0.4% 


0.02% 


0.02% 


0.0017% 



Table 1. Numerical results from ridgerunner compared to the minimum ro- 
pelength values from [ , 11 , | and lITOl . The relative errors in the computations are quite 
small. 




Figure 8. The left-hand picture shows a (loose) configuration of the "simple 
clasp" — a simple two-strand tangle which serves as an interesting model for the 
interaction between two ropes passing over each other at right angles. A ropelength- 
critical configuration of this tangle has been derived and studied extensively in [fTO| 
and [111. Since this derivation included an explicit strut set, it is natural to compare 
ridgerunner's results to this theoretical picture. This comparison is shown in 
the two plots center and right, which plot the positions of struts in arclength coor- 
dinates with the origin located where each curve first begins to turn. The enlarged 
plot (right) shows the agreement between theoretical and computational results. The 
data shown is from a 332 edge polygonal clasp. 



a residual of 0.001 for 202 of the knots and links in our data set, proving that our knots are close 
to being critical for the CThi thickness. While our knots are not quite equilateral, they all satisfy 
the hypotheses of Corollary 2.8 and are hence also close to critical for the original Thip thickness. 
Because of this corollary, we know that both thicknesses are equal for our configurations, so we 
have computed and reported the Thip thickness and ropelength below. 



We started each knot from at least five initial configurations, including the configurations from 
KnotPlot [|49l (similar to the configurations in Rolfsen's table), the TOROS simulated annealer [47], 
Gilbert's minimized configurations from the online Knot Atlas [[24[[. hand-drawn configurations 
from Kawauchi's A Survey of Knot Theory [[29l , and positions generated from KnotPlot's di agr am 
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command. The results shown describe the lowest ropelength we achieved from any of these starting 
configurations. 



The polygonal ropelengths for our curves appear in the column Rop^ of Tables [3| |5|of Appen- 
dix |A[ while a plot of the ropelengths organized by crossing number appears in Figure 10 



Generating upper bounds for smooth ropelength. Our computations yielded a large set of 
approximate minimizers of Lcn(vyThij,(v). From these, we wanted to generate upper bounds on the 
minimum (smooth) ropelength of these knots and links. Rawdon has given general bounds P31l47ll 
on the rate at which Thip — )■ Thi which we could have used for this purpose. But we were 
interested in small improvements in ropelength, so we used a more careful approach. 

Our procedure for constructing smooth ropelength bounds from polygonal data is as follows. 
Beginning with V, we splice circle arcs of radius MinRad(fj) into the comers at vertices Vi as 
shown on the left-hand side of Figure [9] to create a piecewise curve V{s). The minimal radius 
of curvature for this curve is equal to MinRad(V). But the self-distances of V{s) may be different 
from those of the polygon V if they involve the new circle arcs. 



X 





^^+1 / X + (1 + K)e^ 

Figure 9. On the left, we see the curve constructed from splicing a circular arc of 
radius MinRad(t>j) into t>j_it>jt>j_|_i. This curve is C\ but not at the splice points. 



On the right, we see the setup for Proposition 5.1 On the left and right are arcs c 
and d with curvature < K and length < e. The minimum distance x between them 
occurs at c(0), d{0). We prove that the distance between any other pair of points 

c(so) and (i(to) is bounded above by x + (1 + K)e 
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We must therefore compute the self-distances of V{s). This poses a problem: V(s) is composed 
of arcs of circles and line segments and Neff has shown that there is no simple formula for the 
distance between two arbitrary circle arcs in 3-space fi38il . So we estimate the self-distances of 
the smooth curve V{s) by taking distances between a finite number of sample points on the curve 
separated from one another by some e. We bound the error in our computation in terms of e using 
the following Proposition. 

Proposition 5.1. Suppose that c(s) and d{t) are each unit-speed piecewise arcs with curvature 
bounded above by K. Further, suppose that ||c(0) — d{Q) || > 1/2 is the minimum distance between 
c and d. Then for any < Sq, to ^ ^ 

||c(so)-rf(to)|| < ||c(0)-d(0)|| + (l + ir)el 



Proof. Since ||c(s) — d{t) \\ has a local min at (0, 0), we know that 

(c(0),c(0) -rf(0)) = 0, and (d'(0), c(0) - d(0)) = 0. 

25 



Further, the curvature bound tells us that ||c"||, < K. We will use these facts to estimate 

||c(so) - d{tQ)f. If we let C{sq) = J^" c'(s) ds and D{to) = /J" d'{t) dt then we have c(so) = 
C(so) + c(0) and d{to) = D{to) + d{0), so 

(19) ||c(so) - d{to)f = \\Ciso) - Dito)f - 2 (C(so) - Z^(to), c(0) - rf(0)) + ||c(0) - rf(0)f . 

Since c(s) and (i(t) are unit- speed curves, and < SQ,to < e we know that ||C(so)||, 111^(^0)11 < e 
and so the first term is bounded above by 4e^. 

The middle term is more interesting. As before, we can let CC{s) = c"{x) dx and DD{t) = 

Jo d"{y) dy, so c'{s) = CC{s) + c'{0) and d'{t) = DD{t) + d'{0). Since c'(0) and d'{0) are normal 
to c(0) — d{0), we can then write this middle term as 

- (C(so) - D{to), c(0) - d{0)) = -{ CCis) ds - / DD{t) dt, c(0) - d{0) 

\Jo Jo I 

Since ||c"||, \<X'\ < K, we know < Ks, \\DD{t)\\ < Kt. Thus (remembering that sq, 

to < e) the norms of the integrals on the right above are each bounded above by and the 

entire dot product is bounded above by Ke^\\c{Q) — (i(0) ||. 



Thus the right hand side of ^ is bounded by ||c(0) - d(0) f + 46^ + 2Ke'^\\c{Q) - d(0) || . Since 
1/2 < ||c(0) -d(0)||,4e2 < 2e2||c(0) -rf(0)||. Using this, we see that 

4e2 + 2ife2||c(0) - rf(0)|| + ||c(0) - rf(0)f < ||c(0) - d(0)f + (2 + 2K)e^\\c{<d) - ci(0)|| 
< ||c(0) - rf(0)f + (2 + 2ir)e2||c(0) - t/(0)|| + (1 + K)V 

= (||c(0)-rf(0)|| + (l + ir)e^)^ 
This completes the proof. □ 

Our code, named roundout_rl|^ establishes a coarse net of points on V{s) x V{s) ~ [0, 1] x 



[0, 1] and then eliminates subsquares of this square from consideration using Proposition 5.1 The 
remaining squares are then subdivided and searched in turn. The process terminates once we have 
computed the local minima of d{p, q) on the square with whatever accuracy we require. 

Using roundout_rl in double-precision machine arithmetic we found upper bounds for the 
ropelengths of our 379 minimized configurations. These figures appear in column Rop of Tables[3]- 
|5]of Appendix [A} These figures constitute the best known dataset on the lengths of tight knots and 
links. The data is summarized in Figure 10 and Table [2j 



To test how accurate these final results are likely to be, we computed the relative residual 
II 11/11 _v/|| for all these knots and links. The average residual of knots in our tabulation is 
about 0.00299. We have achieved residuals as low as 2.54 x 10^^ for knots and links of special 
interest, such as 818, IO123, the trefoil, and the Borromean rings. A table of these residuals appears 
in Appendix |Aj Four knots and links in our calculation turned out to be particularly difficult for 
ridgerunner: lOei, Sfg, 8| and 



'freely available as part of the octrope library 
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Figure 10. This graph shows the relationship between ropelength (y-axis) and 
crossing number (x-axis) for knots and links in our data set. The bottom lines show 
the bound of Denne et al. [ITTIl for ropelength of a nontrivial knot (horizontal line) 
and Diao's bound [|T8l for ropelength in terms of crossing number (curve). We 
can see that there is a substantial overlap of ropelength values between different 
crossing numbers. This is reflected in Tables [6[j7] of Appendix [A[ which show the 
knots in ropelength order. Table [2] shows the links of least and most ropelength for 
each crossing number. 



Cr 


Rop 


Links 


3 


32.74 


3i 


4 


[40.0122,42.0887] 


4?,4i 


5 


[47.2016,49.7716] 


5l, 5;^ 


6 


[50.5539,58.1013] 


^3' "3 


7 


[55.5095,66.3147] 


7?, 7^ 


8 


[60.5754, 75.2592] 


81 8i 


9 


[66.0311,83.6092] 


q2 q2 


10 


[71.0739,92.3565] 


10i24, 10l23 



Table 2. This table shows the links of smallest and largest minimum ropelength 
for each crossing number (according to our data). Recall that we did not mini- 
mize ten-crossing links, so it is likely that some ten-crossing link has more or less 
ropelength than the IO123 and IO124 knots. 

Generation of tightening animations, pictures, and strut sets. We have saved the minimiza- 
tion runs for each of these knots and links as an animation showing the tightening knot. These 
animations are posted on the web at http : / / www . jasoncantarella . com/mo vs/j 

We have also generated images of the polygonal strut sets and approximately tight configurations 
for each of the 379 knots and links in our data set. Space considerations prevent us from including 
all of this data in this paper, so they are enclosed in the associated Atlas of Tight Links Q . Figure 1 1 
shows a typical page from the Atlas. All of our tight knot and link data, including coordinates for 
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Link Rop 



Rop 



Verts Stmts 



Link Rop 



Rop 



Verts Struts 



7i 



61.4234 61.4067 512 



788 



4f 40.0203 40.0122 400 610 



Figure 1 1 . This figure shows simplified versions of two pages from the Atlas of 
Tight Links for the knot 7i and the link Af. On each page, the top left pictures 
show three views of the link. The triangular graphic shows the struts of the link as 
found by ridgerunner plotted as points (s, t) in arclength coordinates along the 
link. The graph on the bottom of the page shows the curvature of the curve. The 
background of each plot changes color to indicate the change from one component 
to the next. The key along the left-to-right diagonal is given in ropelength units and 
color-coded with the pictures at upper left to show which component is referred to 
by the plot. 



the tight configurations, is publicly available with the publication of this paper. We note that for 
technical reasons, our minimized configurations have thickness close to 1/2 (rather than 1, as in 
the discussion above), and hence their maximum curvature is 2. 



Discovery of symmetric tight knots. An interesting feature of the ropelength function is that 
minimizing ropelength seems to break any symmetry enjoyed by the original configuration of a 
given knot. For instance, while the minimizing configuration for the (3, 2) torus knot 3i appears 
to be threefold symmetric (as expected), the minimizing configuration for the (5, 2) torus knot 5i 
is not fivefold symmetric. It was therefore somewhat surprising to discover two knots in our data 
set, 8i8 and IO123 for which the tight configurations are highly symmetric. These knots are shown 
in Figure [12} Their self-contact sets (which appear on pages 67 and 358 of the Atlas, and are 



reproduced in the Appendix of this paper on pages 36 and 37 1 are highly suggestive, resembling 
those of the Borromean rings (page 29), and appearing to consist of a single element repeated 
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Figure 12. Two highly symmetric tight knots are the Sig knot shown above left 
and the IO123 knot shown above right. Rounding the corners of these curves yields 
ropelength upper bounds of 74.9063 and 92.3565, respectively. Because their strut 
sets break into a particularly simple form (see pages 36 and 37 ), these knots may be 
better candidates for an explicit solution than seemingly simpler knots such as the 
trefoil. 



several times. This feature implies that these knots may be better candidates for explicit solution 
than the seemingly simpler trefoil knot. 



6. Future directions 



Several directions for future research suggest themselves from these experiments. First, we note 
that while we have given finite strut sets for several polygonal knots and observed that they are 
close to the the 1 -dimensional strut sets for the corresponding smooth tight configurations, we 
have not proved a theorem explaining how our polygonal strut sets converge to the strut sets of 
a critical polygon. We conjecture that this is part of a larger theorem which would show that if 
a family of polygonal ropelength critical configurations V„ converge to a C^'^ curve V then V is 
ropelength critical in the sense of [9J, the strut sets of the V„ converge in Hausdorff distance to 
the self-contact set of V, and the kink sets of the Vn converge to the portion of V at maximum 
curvature. 

There are several features of the tight knot data set that we have discovered that seem worthy 
of further investigation. Carlen, Smutny and Maddocks noted in [14J that curvature constraints 
seemed to be "within a rather small tolerance of being active" at several points on their numerical 
approximations of the tight trefoil and figure-eight knots. Baranska et al. provided numerically 
smoothed plots of the curvature of their approximately tight trefoil in ^ which appear to confirm 
this observation. 
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6.55 
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Figure 13. This plot shows the computed 1/MinRad values as a function of 
arclength along the polygon for a 2400 edge trefoil with thickness close to 1/2, 
residual 0.0018 and polygonal ropelength 32.743663 (rounding out the corners as 
described above gives a smooth ropelength upper bound of 32.74352 for this con- 
figuration). The value at each vertex is plotted above with no numerical smoothing. 
Though there is some noise in the portions of the plot where curvature is not con- 
strained, the six kinked regions are clearly resolved. A total of 117 vertices are 
involved in these regions. 



We have noticed the same phenomenon in our data sets. Our computation of the curvature for 
the trefoil appears in Figure [T3j In the Atlas of Tight Knots, we highlight the active curvature 
constraints found by ridge runner as part of the minimization process by red lines on the plot 
of strut sets. These occur in 372 of the 379 knots and links minimized. This provides suggestive 
numerical evidence that kinks are rather common in tight knots. We intend to provide better 
evidence for this conjecture in an upcoming publication. 



Several authors have proved versions of the theorem that an interval of a tight knot with curvature 
less than the maximum allowed and no struts must be a straight line segment ll27ll50l l9ll2T1l. We 
see this phenomenon 338 times in the Atlas, for instance in the link 6| on page 28 of the Atlas (see 
also Figure 14), which appears to have three straight segments of length 2.1, 1.14, and 0.56. We 
highlight these segments in blue on the plots in the Atlas. These segments are almost as common 
as kinked regions in our data set, suggesting that they are generic features of tight configurations. 
Gonzalez has conjectured that every composite knot formed from joining a knot to its mirror image 
has a critical configuration with a pair of straight segments. We do not address this conjecture 
here since we only consider prime knots and links, but we do intend to compute approximately 
minimizing composite knots and links in a future publication. 



The paper (9) (as well as ll34l under very different hypotheses) shows that a pair of arcs in a 
tight knot coparametrized by a single family of struts and having curvature less than the maximum 
bound form a standard double helix. As far as we can tell, this phenomenon only occurs a few 
times in the Atlas, for instance in the 63 link on page 28, the 7y link on page 43, the 819 knot on 
page 66, and possibly in the 87 link on page 91. It would be interesting to look for more critical 
configurations with double -helix sections. 
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Figure 14. This figure shows two views of our computed tight configuration of 
the link 63 (ropelength upper bound 50.5539). Straight segments on the blue and 
white components, which occur when these components lose contact with the other 
components of the link, are highlighted in darker blue. 

We also contemplate further improvements to our numerical knot tightening methods. The con- 
strained gradient descent method presented in this paper is a significant improvement over simu- 
lated annealing — in practice, it has proved to be an effective minimizer for both knots and links. 
But this is surely not the last word in numerical ropelength minimization. Our method is a mem- 
ber of the class of "projected-gradient" methods introduced by Rosen and Zoutendijk in the early 
1960's [|48ll54l . These algorithms are subject to a number of well-known numerical problems, such 
as a tendency to "wobble" when confronted with a steep-sided valley and the problem of "zigzag- 
ging", which occurs when elements repeatedly enter and leave the strut and kink sets on successive 
minimization or error-correction steps. Our implementation seems to suffer from both these prob- 
lems during some difficult minimizations. We have experimented with adding conjugate-gradient 
features to our existing code to solve these problems, but so far the results seem to yield only a 
slight improvement. 

For these reasons, more modern methods such as sequential quadratic programming (SQP) have 
become the norm [|23l . Codes implementing these methods require the user to specify a set of 
constraint functions in advance. Unfortunately, in our formulation of the constraint thickness an n- 
vertex polygon has O(n^) self-distance constraints and 0{n) turning angle or MinRad constraints. 
For a typical polygon with 10"^ vertices, this would mean a set of 10^ constraints — too many to 
be practical. However, if we know approximately which self-distance constraints will be active 
in the final configuration, we can ignore constraints that we expect to be inactive, resulting in a 
reduced constraint set of size 0{n). Our approximately minimized polygons provide exactly this 
information. For this reason we imagine an important use of our data will be in formulating input 
problems for a future SQP-based knot-minimizer. Our polygons are already serving as input for 
the biarc-based annealer of Carlen, Smutny, and Maddocks [fT4l . 

While our data set is detailed and suggestive, solving explicitly for the structure of ropelength 
minimizing (smooth) knots and links is likely to require even better data. Cantarella et al. [9] have 
shown that a critical shape for the simple clasp formed when ropes pass over one another at right 
angles contains tiny straight segments of length a few thousandths of the total length of the curves. 
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Resolving these features will require converged runs for polygonal ropelength minimizers with tens 
of thousands of vertices, an ambitious goal that will keep this area of experimental mathematics 
active for some time to come. 
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Appendix A. Ropelength Data 



The pages that follow contain three sets of tables of ropelength data. The first set, Tables |3|(5] on 



pages 38 -40, show the polygonal ropelength (Ropp) and ropelength upper bounds (Rop) that we 
have obtained for each of the knot types that we have considered. The knots and links are organized 
according to their position in Rolfsen's table, with the link being the z-th example of a prime 
X-crossing link of y components in the table. We have identified the two "Perko pair" knots lOiei 
and 10i62 and renumbered the subsequent knots accordingly, so there are only 165 ten-crossing 
knots in our results. 



The second set. Tables |6] and |7] on pages 41-42 show the same knot and link types ordered by 



ropelength upper bound. These tables are to be read down each column from the top left to the 
bottom right. We can see that this order is quite different from the one in Rolfsen's table with (for 
instance) the 2-component link 7^ occurring before any 6 or 7 crossing knot and the IO124 knot 
occurring before many 8 and 9 crossing links. 



The third set of tables. Tables [8- 10 on pages 43 -45 give the residual of each of our computed 



configurations. The low residuals show that they are close to critical in the sense of Theorem 2.18 
We include this data as measure of the relative quality of each of our minimized configurations. 



On pages 36 -37 are reproductions of the pages from the Atlas of Tight Knots for the approxi- 
mately tight 818 and IO123 knots . On the top left of each page are three views of the tight configu- 
rations, with kinked regions highlighted in red. On the top right is a plot of the self-contact map of 
the configuration. Each of these plots consists of a triangular region with the hypotenuse labeled 
with arclength values on the knot. A green box is plotted at (s, t) on the plot if there is a strut con- 
necting L(s) and L(t). Below the graph appears a plot of 1/ MinRad for the polygon (to the same 
scale). Kinked regions of maximum curvature are plotted in red on the graph. Each such region has 
a key on the right-hand side of the plot showing the arclength positions of the start and end of the 
kink (in order to give a sense of the relative scale of the kinked region). At the bottom of the page 
is a line of data giving the polygonal ropelength Rop^ (as measured by octrope), ropelength 
upper bound Rop (from roundout_rl), filename, number of vertices and struts, maximum and 
minimum curvature values and number of kinked regions. The last entry shows the total arclength 
of straight regions in the curves (0 for these two knots, but nonzero for many knots and links in 
the Atlas). 
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Table 3. Part 1 of Ropelengths of Tight Knots and Links by Knot Type 
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Table 4. Part 2 of Ropelengths of Tight Knots and Links by Knot Type 
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86.8514 


9?7 


80, 


.3223 


80, 


.3022 












1021 


87.0497 


87.0343 












9?8 


81, 


.7563 


81, 


,7461 


9? 


81, 


.1522 


81, 


,1333 


1022 


87.2417 
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88.1319 


88.1153 
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Table 5. Part 3 of Ropelengths of Tight Knots and Links by Knot Type 



Link 


Rop„ 


Rop 


Link 


Ropp 


Rop 


Link 


Ropp 


Rop 


IO37 


88.1319 


88.1153 


lOso 


89, 


.1669 


89, 


,1556 


IO123 


92, 


,3646 


92, 


,3565 


IO38 


88.3478 


88.3257 


1081 


90, 


.0181 


90, 


.0007 


IO124 


71, 


.0894 


71, 


.0739 


IO39 


88.3562 


88.3323 


1082 


88, 


.7011 


88, 


.6801 


IO125 


74, 


.9907 


74, 


.9778 


IO40 


89.2659 


89.2464 


1083 


89, 


.5544 


89, 


,5314 


IO126 


77, 


,6202 


77, 


,6026 


IO41 


89.0725 


89.0553 


1084 


89, 


.6518 


89, 


.6788 


IO127 


80, 


.0235 


80, 


.0124 


IO42 


89.9013 


89.8857 


1085 


87, 


.8403 


87, 


.8164 


IO128 


76, 


.4187 


76, 


.4026 


IO43 


89.3512 


89.3366 


1086 


88, 


.7050 


88, 


,6851 


IO129 


78, 


,5739 


78, 


,5553 


IO44 


88.8714 


88.8515 


1087 


89, 


.1363 


89, 


.1173 


IO130 


78, 


.8499 


78, 


.8356 


IO45 


89.4836 


89.4621 


1088 


89, 


.5638 


89, 


.5461 


IO131 


81, 


.2871 


81, 


.2678 


1046 


86.4718 


86.4487 


1089 


89, 


.4343 


89, 


,4178 


IO132 


74, 


,7441 


74, 


.7330 


IO47 


87.3043 


87.2821 


IO90 


88, 


.9330 


88, 


.9115 


IO133 


77, 


.1813 


77, 


.1631 


IO48 


87.3814 


87.3643 


IO9I 


88, 


.9611 


88, 


.9435 


IO134 


78, 


.6521 


78, 


.6377 


IO49 


88.2914 


88.2705 


IO92 


89, 


.6200 


89, 


,6011 


IO135 


81, 


,2305 


81, 


,2157 


IO5O 


87.3876 


87.3716 


IO93 


88, 


.3962 


88, 


.3773 


10l36 


78, 


.0398 


78, 


.0276 


IO5I 


88.3209 


88.3002 


IO94 


88, 


.8514 


88, 


.8306 


IO137 


79, 


.6352 


79, 


.6185 


IO52 


88.0719 


88.0565 


IO95 


90, 


.0056 


89, 


,9848 


10l38 


82, 


,5504 


82, 


,5320 


IO53 


88.8361 


88.8180 


IO96 


89, 


.5493 


89, 


.5284 


IO139 


72, 


.9001 


72, 


.8944 


IO54 


87.5336 


87.5127 


IO97 


89, 


.4340 


89, 


.4163 


10l40 


73, 


.8610 


73, 


.8477 


IO55 


88.3760 


88.3699 


IO98 


89, 


.7172 


89, 


,6969 


IO141 


76, 


,9687 


76, 


.9543 


IO56 


89.0160 


88.9973 


IO99 


88, 


.8926 


88, 


.8734 


10l42 


75, 


.8951 


75, 


.8754 


IO57 


89.6126 


89.5946 


lOloo 


88, 


.7124 


88, 


,6927 


IO143 


78, 


,2422 


78, 


.2307 


IO58 


88.9623 


88.9445 


lOioi 


89, 


.7344 


89, 


,7210 


IO144 


81, 


,4378 


81, 


,4275 


IO59 


89.2228 


89.2090 


IO102 


88, 


.7969 


88, 


.7734 


IO145 


75, 


.9194 


75, 


.9076 


lOeo 


89.3397 


89.3190 


IO103 


88, 


.7971 


88, 


.7914 


10l46 


79, 


.7416 


79, 


.7322 


lOei 


86.4755 


86.4561 


IO104 


91, 


.7476 


91, 


,7280 


IO147 


79, 


,1666 


79, 


,1571 


1062 


87.5318 


87.5071 


IO105 


89, 


.8260 


89, 


.8055 


10l48 


79, 


.0893 


79, 


.0742 


1063 


88.4046 


88.3861 


IO106 


89, 


.1546 


89, 


,1319 


IO149 


81, 


.0500 


81, 


.0318 


1064 


87.4878 


87.4742 


IO107 


89, 


.7525 


89, 


,7356 


IO15O 


80, 


,1392 


80, 


,1219 


1065 


88.3918 


88.3725 


IO108 


88, 


.5137 


88, 


.4932 


IO151 


81, 


.8414 


81, 


.8207 


lOea 


89.0275 


89.0047 


IO109 


91, 


.1966 


91, 


.1789 


IO152 


79, 


.1715 


79, 


.1556 


1067 


88.4741 


88.4534 


lOiio 


89, 


.6275 


89, 


,6114 


IO153 


80, 


,4764 


80, 


,4648 


1068 


88.1199 


88.1013 


lOiii 


89, 


.6677 


89, 


.6438 


IO154 


81, 


.5405 


81, 


.5218 


1069 


89.0983 


89.0778 


IO112 


89, 


.5744 


89, 


.5529 


IO155 


78, 


,0648 


78, 


.0503 


IO70 


89.2068 


89.1846 


IO113 


90, 


.2239 


90, 


,2141 


10l56 


79, 


,5639 


79, 


,5443 


IO71 


89.0853 


89.0699 


IO114 


89, 


.3062 


89, 


.2856 


IO157 


81, 


.4731 


81, 


.4568 


If) -70 


89.1974 


89.1779 




90, 


.4340 


90, 


.4176 


-*-'-'158 


81, 


.6398 


81, 


.6220 


IO73 


89.5332 


89.5130 


IO1I6 


90, 


.2703 


90, 


,2583 


IO159 


79, 


,8863 


79, 


.8692 


IO74 


88.1285 


88.1077 


IO117 


89, 


.5335 


89, 


.5245 


IO16O 


78, 


.1529 


78, 


.1472 


IO75 


88.9725 


88.9524 


IO1I8 


89, 


.5261 


89, 


.5094 


IO16I 


74, 


.5460 


74, 


.5302 


1076 


88.3673 


88.3479 


IO119 


90, 


.1394 


90, 


,1226 


10l62 


81, 


,0033 


80, 


,9838 


IO77 


88.5689 


88.5471 


IO12O 


90, 


.1862 


90, 


.1674 


10l63 


82, 


.6629 


82, 


.6548 


IO78 


88.5548 


88.5322 


IO12I 


89, 


.9375 


89, 


.9240 


10l64 


82, 


.1862 


82, 


.1698 


IO79 


88.9647 


88.9488 


IO122 


89, 


.8258 


89, 


,8094 


10l65 


82, 


.8211 


82, 


.8040 


lOso 


89.1669 


89.1556 


IO123 


92, 


.3646 


92, 


.3565 
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Table 6. Part 1 of Knot and Link Types sorted by Ropelength 



Link 


Link 


Link 


Link 


Link 


Link 


Link 


Link 


2\ 


81 


948 


93 


915 


9l8 


IO21 


IO23 


3i 


o3 
^10 


o2 
^10 


IO143 


r»2 

931 


r\2 

955 


IO12 


lOioo 




946 


q2 
814 


94 


n2 
921 


r>2 

924 


lOio 


IO25 


4i 


%o 


8i 


IO129 


/^9 

936 


938 


IO22 


IO102 


5i 


r\2 

961 


814 


9i 


9l4 


IO151 


IO7 


IO103 


52 


942 


o2 
012 


10134 


9l5 


939 


IO47 


IO53 


r 2 


n2 

9I7 


955 


95 


9l9 


r\2 

934 


IO48 


IO94 


63 


951 


n2 
9^6 


r>2 

9| 


926 


r>2 
939 


IO50 


IO44 


61 


81 


10l40 


IO130 


9i 


937 


IO15 


IO99 


7? 


n2 

9g4 


815 


10l48 


924 


r\2 

9l2 


IO16 


IO13 


6i 


82 


949 


IO152 


91 


9lo 


1064 


IO90 


61 


9I5 


o2 

85 


IO147 


IO162 


9ii 


1028 


IO91 


62 


1O124 


812 


935 


r\2 

9^7 


r\2 

9?i 


1062 


IO58 


^2 


83 


948 


n2 

9j3 


IO149 


9? 


IO19 


IO79 


61 


82 


945 


9l9 


922 


926 


IO54 


IO75 


63 


84 


9l7 


r\9 

9i 


9^2 


929 


1085 


IO56 


6i 


943 


o2 
Oi3 


10l56 


921 


933 


IO36 


1066 


63 


944 


rv2 

958 


96 


91 


10l64 


IO52 


IO41 


87 


85 


n3 
9f5 


10137 


925 


934 


IO35 


IO71 


819 


r\2 

944 


o2 

87 


r\2 

95 


9? 


r\2 

930 


1068 


1069 


7i 


957 


9l4 


97 


9? 


9§ 


IO74 


1087 


820 


9?3 


817 


10l46 


928 


9?o 


IO37 


IO106 


72 


87 


IO161 


9io 


937 


r>3 
9l2 


IO3I 


1080 


73 


81 


85 


10159 


927 


9| 


IO49 


IO72 


^2 

7f 


n2 

9i9 


10132 


99 


929 


10l38 


IO33 


IO70 


74 


9l8 


81 


93 


IO135 


10l63 


IO51 


IO59 


815 


89 


947 


9l4 


r»9 

935 


938 


IO14 


IO29 


8i 


86 


9ii 


98 


IO131 


933 


IO38 


IO40 


1-7? 

7i 


o2 

Si 


818 


10127 


923 


10l65 


IO39 


IO114 


74 


9l9 


10125 


9l2 


928 


941 


IO76 


lOeo 


7.5 


o2 

8i 


84 


n2 

9i 


941 


r\2 

942 


IO3O 


IO43 


7i 


83 


9l6 


9l6 


9l6 


lOi 


IO55 


IO97 


821 


813 


81 


r\2 

920 


932 


IO3 


1065 


1089 


77 


9i2 


9i 


9ii 


IO144 


102 


IO93 


IO45 


76 


^56 


IO142 


IO150 


IO157 


104 


1063 


IO118 


7i 


IO139 


IO145 


936 


930 


108 


IO24 


IO73 


949 


811 


920 


920 


932 


IO46 


IO18 


IO117 


'5 


°2 


J- '-'128 


^23 


1 Di c/ 
J-U154 


J-^61 


1 Ona 


-■-"96 


9I3 


810 


10141 


9i 


9i 


IO5 


1067 


1083 


8! 


8?i 


10133 


9l3 


9l8 


lOe 


IO1O8 


1088 


7i 


8i 


10126 


^9 


9f 


IO9 


IO78 


IO112 


1 


8i 


10l36 


q2 

yi7 


10l58 


IO20 


IO77 


IO57 


816 


8i 


IO155 


q2 

^10 


940 


lOii 


IO32 


IO92 


8^ 


816 


92 


9l7 


931 


IO34 


1082 


lOiio 


9^3 


^60 


IO16O 


IO153 


9i 


IO17 


1086 


lOiii 


8? 


^48 


93 


q2 

^15 


9?8 


IO21 


IO23 


1084 
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Table 7. Part 2 of Knot and Link Types sorted by Ropelength 



Link 


Link 


Link 


Link 


Link 


Link 


Link 


1084 


10l07 


IO27 


IO95 


IO120 


IO115 


IO123 


1098 


10l05 


IO42 


lOsi 


IO113 


IO109 




lOlOl 


10l22 


IO121 


IO119 


IO116 


IO104 




10l07 


1027 


IO95 


IO120 


IO115 


IO123 
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Table 8. Part 1 of Residuals of Tight Knots and Links by Knot Type 



Link 


Residual 


92 
^1 


4 CI 94„ rii^ 


3i 


0.00621792 


4i 


0.000996335 




U.UUUyyyo4y 


5i 


0.00981995 


r 

02 


u.uuyy4/ /o 


5? 


0.00998078 


61 


0.000999592 


62 


0.00897204 


63 


0.000979541 


6i 


0.000999952 


c2 
62 


0.000999833 


c2 
63 


c\ r\c\{\c\{\c\c\ A 

0.00999004 


6? 


0.00998537 


62 


0.000705159 


6i 


0.00627026 




0.00105833 


72 


0.00998149 


73 


0.00999358 


74 


0.00100877 


75 


0.000999532 


76 


0.000979869 


77 


0.00100393 


7? 


0.000999487 


71 


0.00101952 


'3 


0.000999871 


-72 

7| 


0.00099954 


7^ 


0.000999894 


72 

'6 


0.00100556 


72 
'7 


0.00320787 


71 


0.0018494 


7? 


0.000999748 


81 


0.00898769 


82 


0.000982684 



Link Residual 



Link Residual 



Link Residual 



82 


0.000982684 


83 


0.00100028 


84 


0.00100103 


85 


0.00100033 


86 


0.000999848 


87 


0.00101551 


88 


0.000981272 


89 


0.000999932 


810 


0.000978418 


811 


0.000979921 


812 


0.00998976 


813 


0.000993117 


814 


0.000981486 


815 


0.0099948 


816 


0.000981316 


817 


0.00999085 


818 


0.000900015 


819 


0.000998339 


820 


0.00099998 


821 


0.000999988 


8? 


0.00100142 


81 


0.000979836 


Si 




0.000999961 


81 


0.00216462 


85 


0.00999516 


8i 


0.00100295 


8? 


0.000999802 


81 




0.000999762 


8i 


0.000979774 


810 


0.000999858 


8fi 


0.00997927 


812 


0.000999968 


8?q 


0.0010008 


814 


0.00101123 


815 


0.00099994 


816 


0.000997563 


8? 


0.00100589 


81 


0.000999904 


8i 


0.00100014 


81 


0.00999606 


8i 


0.000995844 


8i 


0.00099824 


83 


0.00119532 


8i 


0.00100655 



8i 


0.00100655 


si 


0.000980533 


810 


0.0208108 


s\ 


0.00100006 


s\ 


0.000999682 


s\ 


0.780186 


9i 


0.00802077 


92 


0.00997484 


93 


0.00998254 


94 


8.64059e - 05 


95 


0.00999417 


96 


0.000980197 


97 


0.000979897 


98 


0.00101007 


99 


0.000999938 


9io 


0.00113523 


9ii 


0.000981742 


9l2 


0.000979842 


9l3 


0.00999582 


9l4 


0.000984327 


9l5 


0.000979831 


9l6 


0.000999818 


9i7 


0.00100032 


9l8 


0.00992217 


9l9 


0.000981217 


920 


0.00100005 


921 


0.0010001 


922 


0.000998846 


923 


0.000979562 


924 


0.000999907 


925 


0.000977105 


926 


0.00100048 


927 


0.00999324 


928 


0.00996501 


929 


0.000979844 


930 


0.000979942 


931 


0.000979062 


932 


0.000997746 


933 


0.00100114 


934 


0.000999697 


935 


0.000981383 


936 


0.000978472 


937 


0.00999228 


938 


0.000978978 



^38 


0, 


.000978978 


939 


0, 


.000999482 


940 


0, 


.000999343 


941 


0, 


.00899161 


942 


0, 


.000999996 


943 


0, 


.00898749 


944 


0, 


.000999789 


945 


0, 


.0099754 


946 


0, 


.00099973 


947 


0, 


.000998991 


948 


0, 


.00998933 


949 


0, 


.00099957 


9? 

1 


0, 


.00107787 


9i 


0, 


.00100115 


92 




0, 


.00100055 


9l 


0, 


.00099991 


9i 


0, 


.00100118 


9i 


0, 


.00126944 


9? 


0, 


.00104121 




0, 


.00100133 


91 


0, 


.000999724 


92^ 

10 


0, 


.00140283 


9?i 


0, 


.000999221 


9?2 
iz 


0, 


.00100137 


9?3 
lo 


0, 


.00100112 


9?4 

14 


0, 


.000999788 


9?5 
lb 


0, 


.000999236 


9?fi 

lo 


0, 


.00605 


9?7 

1 / 


0, 


.00899775 


9?8 

" lo 


0, 


.000999648 


9?9 
ly 


0, 


.00100405 


9in 


0, 


.000999853 




0, 


.00898977 


9I2 


0, 


.00943088 


913 


0, 


.000998181 


914 


0, 


.000999946 


q2 

^25 


0, 


.0009999 


q2 

^26 


0, 


.00100243 


q2 

^27 


0, 


.00099997 


q2 

^28 


0, 


.000998883 


q2 

^29 


0, 


.00100157 


q2 

^30 


0, 


.00099989 


q2 

^31 


0, 


.000999523 


q2 

^32 


0, 


.00100012 


q2 

^33 


0, 


.000999711 
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Table 9. Part 2 of Residuals of Tight Knots and Links by Knot Type 



Link Residual Link Residual 



q2 

33 


0, 


.000999711 


9l6 


0.000999575 


q2 

34 


0, 


.00100169 


9?7 


0.247874 


q2 

35 


0, 


.000999778 


9l8 


0.000999841 


Sie 


0, 


.00100172 


q3 


0.00101035 


Sir 


0, 


.000999058 


920 


0.00100002 


q2 


0, 


.000999748 


921 


0.00100039 


q2 

39 


0, 


.000999888 






9lo 


0, 


.000999835 


9t 


0.000979958 


9li 


0, 


.00100037 






q2 


0, 


.000998679 


lOi 


0.00101691 


q2 

^43 


0, 


.00100109 


IO2 


0.00100023 


q2 


0, 


.00100838 


IO3 


0.000991435 


q2 


0, 


.00997492 


IO4 


0.00100846 


9l6 


0, 


.00100042 


10.; 


0.00100194 


9l7 


0, 


.00999831 


lOe 


0.000979506 


Qls 


0, 


.000999984 


IO7 


0.0097283 


9l9 


0, 


.000999984 


lOg 


0.000980356 


q2 

50 


0, 


.000999226 


IO9 


0.000979784 


q2 


0, 


.000999443 


lOio 


0.00999688 


q2 


0, 


.000999958 


lOii 


0.00760935 


9i3 


0, 


.00996962 


1012 


0.000991292 




0, 


.000999703 


1013 


0.000999947 


q2 

55 


0, 


.00100064 


1014 


0.0010261 


9i6 


0, 


.000979788 


1015 


0.000979185 




0, 


.00255237 


1016 


0.000985699 


q2 

^58 


0, 


.000999155 


1017 


0.00998848 


q2 

^59 


0, 


.00108631 


1018 


0.000979621 


q2 

60 


0, 


.000999312 


1019 


0.00098045 


q2 

61 


0, 


.00100091 


1020 


0.000979959 








1021 


0.000999057 


9? 


0, 


.000999763 


1022 


0.000991413 


9i 


0, 


.000999746 


1023 


0.00999682 


9i 


0, 


.00100525 


1024 


0.00166886 


4 


0, 


.000999641 


1025 


0.000994731 


9i 


0, 


.00100042 


1026 


0.00098015 


9i 


0, 


.000999746 


1027 


0.000999869 


9? 


0, 


.000999935 


1028 


0.00996703 


9i 


0, 


.000999751 


1029 


0.00116525 


"9 


0, 


.000996684 


1030 


0.000999376 


q3 


0, 


.00099985 


1031 


0.000979897 


q3 
"ii 


0, 


.0010755 


1032 


0.000979993 


q3 

^12 


0, 


.00100439 


1033 


0.000979857 


q3 

"l3 


0, 


.00980919 


1034 


0.00098555 


q3 


0, 


.00900147 


1035 


0.000982115 


q3 

^15 


0, 


.00112426 


1036 


0.000979692 


q3 

^16 


0, 


.000999575 


1037 


0.000999835 



Link 


Residual 


Link 


Residual 


IO37 


0, 


.000999835 


1082 


0, 


,000978946 


IO38 


0, 


.000979821 


1083 


0. 


.00999433 


IO39 


0, 


.000986038 


1084 


0. 


.0099812 


IO40 


0, 


.00100863 


1085 


0, 


,000981325 


IO41 


0, 


.00999693 


1086 


0. 


.000978499 


IO42 


0, 


.000999751 


1087 


0. 


.000979621 


IO43 


0, 


.000980157 


1088 


0, 


,000979845 


IO44 


0, 


.00322255 


1089 


0. 


,0010019 


IO45 


0, 


.000982692 


IO9O 


0. 


,000980234 


1046 


0, 


.00997656 


IO91 


0, 


,000977397 


IO47 


0, 


.000980999 


IO92 


0. 


,00100005 


IO48 


0, 


.00999602 


IO93 


0. 


,000979652 


IO49 


0, 


.000998073 


IO94 


0, 


,00097991 


IO5O 


0, 


.000981787 


IO95 


0. 


,000979668 


IO5I 


0. 


.00098231 


IO96 


0. 


,00018365 


IO52 


0, 


.000999419 


IO97 


0, 


,000999872 


IO53 


0, 


.00101025 


IO98 


0. 


,00999481 


IO54 


0, 


.00999263 


IO99 


0. 


,0099926 


IO55 


0, 


.00998728 


lOloo 


0, 


,00101003 


IO56 


0, 


.00999185 


lOioi 


0. 


,00999705 


IO57 


0. 


.000999798 


IO102 


0. 


,000979674 


IO58 


0, 


.000999966 


IO103 


0, 


,00999479 


IO59 


0, 


.00995441 


IO104 


0. 


,00999683 


lOao 


0. 


.000980266 


IO105 


0. 


,000979902 


lOei 


0, 


.0241498 


IO106 


0, 


,000979055 


1062 


0, 


.00105699 


IO107 


0. 


,000980096 


1063 


0. 


.00998227 


IO108 


0. 


,00127554 


1064 


0, 


.00997603 


IO109 


0, 


,000979798 


1065 


0, 


.00135295 


lOiio 


0. 


,000979638 


1066 


0. 


.000999872 


lOiii 


0. 


,000979851 


1067 


0, 


.000979823 


IO112 


0, 


,00104599 


1068 


0, 


.00100695 


IO113 


0. 


,00999934 


1069 


0, 


.000999786 


IO114 


0. 


,00100087 


IO70 


0, 


.000980057 


IO115 


0, 


,000978725 


IO71 


0, 


.00999226 


IO116 


0. 


,00998661 


IO72 


0, 


.000999942 


IO117 


0. 


,00998396 


IO73 


0, 


.00998888 


IO118 


0, 


,00099987 


IO74 


0, 


.000978382 


IO119 


0. 


,000999834 


IO75 


0. 


.000981812 


IO120 


0. 


,00100037 


1076 


0, 


.000980892 


IO121 


0, 


,00099989 


IO77 


0, 


.00999768 


IO122 


0. 


,000999203 


IO78 


0, 


.000981017 


IO123 


0. 


,0016528 


IO79 


0, 


.0010001 


IO124 


0, 


,00100133 


lOso 


0, 


.000979926 


IO125 


0, 


,00998345 


lOsi 


0, 


.000981576 


IO126 


0. 


,00999723 


1082 


0. 


.000978946 


IO127 


0. 


,00998882 
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Link 


Residual 


Link 


Residual 


Link 


Residual 


Link 


Residual 


IO127 


0.00998882 


IO137 


0.000979856 


IO147 


0.000999813 


IO157 


0.000979535 


IO128 


0.000988223 


10l38 


0.00899453 


10l48 


0.000981385 


10l58 


0.000980822 




00902523 


IO1QQ 
J- '-'139 


000979731 


IO1/IQ 
-*-'-'14y 


00100026 


IO1 KQ 

J- '-'159 


000979791 


IO130 


0.000999987 


10l40 


0.0099924 


IO150 


0.000979903 


IO16O 


0.00998455 


IO131 


0.00959976 


IO141 


0.00100144 


IO151 


0.000979813 


IO16I 


0.00899311 


IO132 


0.000980876 


10l42 


0.000980204 


IO152 


0.00999625 


10l62 


0.000985909 


IO133 


0.000980018 


IO143 


0.00993363 


IO153 


0.0091785 


10l63 


0.00899697 


IO134 


0.00999485 


IO144 


0.00995796 


IO154 


0.00115132 


10l64 


0.000979519 


IO135 


0.00100006 


IO145 


0.00102699 


IO155 


0.00998753 


10l65 


0.000979783 


10l36 


0.00999149 


10l46 


0.00998505 


10l56 


0.0009799 






IO137 


0.000979856 


IO147 


0.000999813 


IO157 


0.000979535 
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